Transcriptome and metabolome analysis revealed the dynamic change of bioactive compounds of Fructus Ligustri Lucidi

Background The Fructus Ligustri Lucidi, the fruit of Ligustrum lucidum, contains a variety of bioactive compounds, such as flavonoids, triterpenoids, and secoiridoids. The proportions of these compounds vary greatly during the different fruit development periods of Fructus Ligustri Lucidi. However, a clear understanding of how the proportions of the compounds and their regulatory biosynthetic mechanisms change across the different fruit development periods of Fructus Ligustri Lucidi is still lacking. Results In this study, metabolite profiling and transcriptome analysis of six fruit development periods (45 DAF, 75 DAF, 112 DAF, 135 DAF, 170 DAF, and 195 DAF) were performed. Seventy compounds were tentatively identified, of which secoiridoids were the most abundant. Eleven identified compounds were quantified by high performance liquid chromatography. A total of 103,058 unigenes were obtained from six periods of Fructus Ligustri Lucidi. Furthermore, candidate genes involved in triterpenoids, phenylethanols, and oleoside-type secoiridoid biosynthesis were identified and analyzed. The in vitro enzyme activities of nine glycosyltransferases involved in salidroside biosynthesis revealed that they can catalyze trysol and hydroxytyrosol to salidroside and hydroxylsalidroside. Conclusions These results provide valuable information to clarify the profile and molecular regulatory mechanisms of metabolite biosynthesis, and also in optimizing the harvest time of this fruit. Supplementary Information The online version contains supplementary material available at 10.1186/s12870-024-05096-3.


Background
Fructus Ligustri Lucidi (FLL) is the fruit of Ligustrum lucidum, which belongs to the Oleaceae family.Ligustrum lucidum is widely distributed in China, South Korea, India, and Australia and is often used for environmental greening.Its berry fruits, which are used as phytomedicine, are generally ovate, elliptic, or kidneyshaped, with a dark purple or grey-black color; they are usually harvested in winter (Fig. S1).In China, dried mature FLL is named Nvzhenzi, and it is used in traditional Chinese medicine (TCM) [1].Modern chemical research has shown that triterpenes, phenylethanoid glycosides, secoiridoids, and flavonoids are the main secondary metabolites of FLL.These metabolites have a wide range of pharmacological effects, and are often used in clinical applications [2,3].FLL is therefore a promising nutritional food additive for humans.However, it is unclear which harvest time is best for optimizing the amount of active components in FLL, and thus its medicinal value.In addition, not enough is known about the biosynthetic pathways involved, or the dynamic changes of active components to choose the harvest time and utilization of FLL as a food additive.
Based on chemical and pharmacological studies of FLL, its main bioactive secoiridoids are oleuropein, nuezhenoside G13, and specneuzhenide [1][2][3][4][5][6].There has been some researches related to secoiridoid biosynthesis in the Oleaceae family, specifically in olives (Olea europaea L.).Oleoside-type secoiridoids, which are found exclusively contained in the Oleaceae family and the genus Caiophora (Loasaceae), are some of the most frequently isolated and characterized biochemical components of FLL [4,7].Although the biosynthetic pathway of secoiridoids has not been elucidated in FLL, it has been examined in olives for oleoside-type secoiridoids such as oleuropein.That pathway can be described as follows.First, 7-deoxyloganic acid is hydroxylated; this hydroxylation is catalyzed by 7-deoxy-loganic acid 7-epihydroxylase (7eDLH), oleoside methyl ester synthase (OMES), and secoxyloganin synthase (SXS) to produce oleoside-11-methyl ester (OME).OME is then glycosylated by glycosyltransferase (GT) to produce 7-β-1-Dglucopyranosyl-11-methyl oleoside (OME-glu).Finally, OME-glu is then conjugated with phenolic components such as tyrosol and salidroside to produce oleoside-type secoiridoids [8].This biosynthetic pathway may share similarities with the corresponding pathway in FLL.However, the biosynthesis of its major compounds in FLL, as well as transcriptional-level variations during FLL periods, have not been reported.
In this study, six fruit development periods of FLL including 45 days after flowering (DAF), 75 DAF, 112 DAF, 135 DAF, 170 DAF, and 195 DAF were selected for metabolic and transcriptomic analyses to explore the dynamic changes in FLL active components and the functions of key genes.A quantitative analysis of the main compounds was performed using high performance liquid chromatography (HPLC).The compounds and genes related to the biosynthesis of secoiridoids, flavonoids and triterpenoids were analyzed by combining metabolites with the transcriptome.Based on the GTs of Rhodiola rosea, we identified nine homologs in FLL and showed that these GTs catalyze the formation of salidroside and hydroxysalidroside from tyrosol and hydroxytyrosol, respectively.The results revealed that changes in their compound contents and in the expression of FLL-related genes complied with certain rules across the different fruit development periods These results provide a basis for the efficient utilization of FLL across different seasonal collection times and diversified development and utilization of FLL.

Metabolite accumulation during FLL development
The metabolites present in six fruit development periods of FLL were analyzed via LC-MS/MS, and the mixture of the extracts from each sampling period was used as a quality control sample to identify the compounds.The total ion chromatogram of the quality control mixture in the negative and positive modes is exhibited in Fig. S2.Based on the standards and previous researches, a total of 70 compounds were identified, of which secoiridoids (40%), flavonoids (14%), and phenylethanoids (11%) were the most abundant.Details of the identified metabolites are listed in Table 1.The MS/MS spectra of compounds identified using the standards are shown in Fig. S3.The PCA analysis showed that 36 samples (six fruit development periods with six replicates) were separated.The samples in each period were grouped into one cluster, showing a trend of separation among metabolic groups (S4A).
To explore the variations in the main compounds across the development periods of FLL, tyrosol, salidroside, hydroxysalidroside, hydroxytyrosol, specnuezhenide, oleuropein, nuezhenoside G13, ligstroside, ligustroflavone, oleoside-11-methyl ester, and 7-β-1-Dglucopyranosyl-11-methyl oleoside were quantified using HPLC based on a previous method [9], with dry FLL being used as a control.Spenuezhenide was the standard of the Chinese Pharmacopoeia (2020), and the abundance of this compound was 4.3% at 112 DAF.Across sampling periods, spenuezhenide accumulated rapidly at 75 DAF (Fig. 1), reaching the highest level at 112 DAF, and then decreasing gradually.Oleuropein accumulated at 45 DAF, and then decreased gradually, after which it increased at 195 DAF.Nuezhenoside G13 accumulated rapidly at 112 DAF, reached its highest level at 135 DAF (2.5343%), then decreased gradually.The dry FLL had a high content of nuezhenoside G13 and specnuezhenide.The oleuropein content in dry FLL was extremely low (0.0272%), while its abundance in fresh FLL was significant (0.8945-3.1503%FLL was lower than that in the six periods.The trends in the LC-MS-based intensities of the compounds mentioned above, across development stages, were similar to their variation in content (Fig. S4B).

Transcriptome analysis and annotations of different FLL development periods
To explore the underlying mechanisms for the changes in the levels of the major compounds at the gene level, transcriptome analysis of six FLL development periods was performed.The analysis generated 112.3 Gb of clean data, and the Q30 base percentage was more than 93% (Table S3).A total of 103,058 unigenes were assembled, with the length of the unigenes ranging from 301 to 16,743 bp.All unigenes were annotated using common databases, including the NR (NCBI non-redundant protein sequences), NT (NCBI nucleotide sequences), KO (euKaryotic Ortholog groups), SwissProt, PFAM, GO (Gene Ontology), and KEGG (Kyoto Encyclopedia of Genes and Genomes), to which approximately 48, 45.45, 15.36, 32.52, 32.12, 32.12, and 8.24 of the unigenes were mapped, respectively (Table S4).Inter-sample correlation analysis of the 18 samples showed good biological duplication within groups (Fig. S5).Differentially expressed genes (DEGs) were selected using the parameters |log2Fold Change| ≥ 1 and FDR < 0.05, with 36,540 unigenes were screened (Table S5).There were 25,783 S6).

Analysis of the triterpenoid biosynthesis pathway
Triterpenoids, such as oleanolic acid, are important secondary metabolites of FLL.Triterpenoids are generated from 2,3-oxidosqualene, which is mainly derived from the cytosolic mevalonic acid (MVA) pathway and, to a smaller extent, the plastidial methylerythritol phosphate (MEP) pathway.In the MVA pathway, the conversions of isopentenyl diphosphate (IPP .IPP and DMAPP are then catalyzed to FPS by farnesyl pyrophosphate synthase (FPPS).Finally, 2,3-oxidosqualene can be obtained using squalene monooxygenase (SQE) and squalene synthase (SQS) catalysis (Fig. 2).The precursor produced can then be consumed by beta-amyrin synthase (BAS) and CYP716A to produce oleanolic acid.The expression levels of these genes (Fig. 2) showed that the genes related to the MEP and the MVA pathways were highly expressed at 45 DAF, 75 DAF, and 112 DAF, when the color of FLL was still green.The expression of downstream genes like FPPS, SQS, SQE, BAS, and CYP716A increased as fruit developed, suggesting that metabolites were gradually catalyzed
TyDC and AAS belong to the AAAD family of proteins.Owing to the high amino acid sequence identity between individual AAAD paralogs, the annotation of TyDC and AAS can be confusing.Genes such as Cluster-13715.45907,Cluster-13715.20262, and Cluster-137 15.28839 were identified as TyDC/AAS.The MAO was not annotated or blasted in our RNA-seq data, while PAR expression was also found to be similar to TyDC/ AAS expression.Based on T4GT and T8GT sequences in R.  4).Cluster-13715.43413 and Cluster-13715.51696 were highly expressed in 45 DAF, and then their expression gradually decreased.Cluster-13715.7146 was highly expressed at 170 DAF (Fig. S8C).It seems that the variation in the contents of salidroside and hydroxylsalidroside may be caused by the differential expression of these GTs over the six months.
We examined the biochemical activity of the above T8GTs and T4GTs above with recombinant enzymes expressed in E. coli, using tyrosol or hydroxyltyrosol as substrates.HPLC analyses of the enzymatic activity assay revealed that the recombinant enzyme 67,869 could completely metabolize tyrosol and hydroxyltyrosol completely to unknown compounds 1 and 2, respectively, with an Rt of 10.82 min and 7.79 min, respectively; 18,726 reacted with tyrosol and hydroxyltyrosol completely to produce salidroside and hydroxylsalidroside; 48,377 catalyzde most of the tyrosol to salidroside, but could not catalyze hydroxytyrosol; 41,400 and 51,596 converted tyrosol to a lesser extent to compounds near icariside D2 with an Rt of 5.75 min (compound 3), while they did not show activity with hydroxyltyrosol; 35,728 and 48,377 converted tyrosol to salidroside and compound 3, respectively, while they could not react with hydroxyltyrosol; 43,413 did not show activity with tyrosol or hydroxyltyrosol; and 7146 converted tyrosol completely to an unknown compound (compound 4) with an Rt of 10.74 min, but could not catalyze hydroxytyrosol (Fig. 5).C4-glycosylation was not observed in the catalyzation of tyrosol or hydroxyltyrosol.The compound 1, 2, 3 and 4 were needed further analysis.Based on these results, we identified three UGTs (18,726, 35,728, and 48,377) with regio-specific T8GT activity.No GTs with region-specific T4GT activity were found in FLL.

Discussion
Dry FLL is known as a phytomedicine and its plant is a greening plant in many countries.FLL is rich in various bioactive compounds, such as triterpenes, phenylethanoid glycosides, secoiridoids, and flavonoids.However, little is known about the dynamic changes of these main compounds, or their biosynthesis pathway as FLL ripens.Furthermore, considering the high price of olive oil [18], the scarcity of Rhodiola plants [19], and the large amount of FLL wasted every year, a clear understanding of the accumulation rules and molecular mechanism of the active components would be helpful in developing and utilizing FLL for use in food, medicine, and other fields.In this study, 70 compounds were identified across the six periods of FLL, and the main compounds, such as specneuzhenide, oleuropein and nuezhenoside G13, were quantified using HPLC.By combining transcriptome and metabolome analyses, changes in the expressions of genes and metabolites related to flavonoids, phenylethanols, triterpenoids, and secoiridoids were analyzed.The functions of C8-specific GTs in secoiridoids and phenylethanol biosynthesis were characterized using enzymatic activity assay.
Terpenoids, mainly secoiridoids and triterpenoids, are the major secondary metabolites of FLL, and they have a range of biological activities.Triterpenoids are enriched in FLL, reaching up to 4-5%, especially oleanolic acid, with contents ranging from 7.37 to 13.3 mg/g [1].Because of this, FLL is the main source of oleanolic acid in the market [2].In our study, the ion intensity of oleanolic acid was stable across the different periods.Oleanolic acid biosynthesis has been demonstrated due to its various pharmacological effects [20].Oleanolic acid can be generated from 2,3-oxidosqualene, which is generated from FPS by SQS and SQE.Then, 2,3-oxidosqualene is processed by BAS and CYP716A to produce oleanolic acid.The genes related to the MEP and MVA pathways were tended to be expressed in the 45 DAF and 75 DAF period.A large number of precursors substance of terpenoids may have been accumulated in these two sampling times.FPPS was highly expressed from 45 DAF to 135 DAF, guiding the compounds flow to triterpenes.These downstream genes, such as SQE, SQS, BAS, and CYP716A, were expressed in all the stages, and may be the reason for the stable oleanolic acid content, which was consistent with previous reports [21].
Secoiridoids, such as those found in FLL and olives, are abundant in Oleaceae and are widely used in food, lumber, cosmetics, edible olive oil, and medicine [22].Oleuropein, the most abundant compound in olive oil, was the first oleoside-type secoiridoids reported [10].In FLL, elevated levels of oleuropein were observed during the immature fruit periods, specifically 45 DAF, 75 DAF, and 112 DAF, but decreased as the fruit underwent physiological development.These observations align with what has been observed in olives.The extent to which oleuropein degradation occurs may be influenced by the activity of β-glucosidase [4].β-glucosidase activity may have been higher at 135 DAF and 170 DAF than in the immature fruit stages, so the oleuropein content in these two stages was lower.However, the lowest oleuropein content was observed in dry FLL.It is possible that the transition from fresh fruit to dry FLL may affect the stability of this compound.Oleuropein was also enriched in olive oil that was extra-virgin cold-pressed rather than extracted, probably because it is heat-unstable.Due to their similar structures, ligstroside was also likely to undergo degradation by β-glucosidase, resulting a significant accumulation of its content at 45 DAF, 75 DAF, and 112 DAF, followed by a substantial decrease in the subsequent 135 DAF, 170 DAF and 195 DAF.In addition, specnuezhenide and nuezhenoside G13 both share a common β-glucose structure and are susceptible to the activity of β-glucosidase.This enzymatic action may lead to a reduction in their contents at 170 DAF and 195 DAF.In contrast, these two compounds exhibited higher levels of enrichment in dry FLL (fully ripened and dried fruits) compared to fresh fruits at 170 DAF and 195 DAF.It is possible that the activity of β-glucosidase was inhibited during the drying process, resulting in the cessation of the conversion of these two products.Additionally, the activity of catalytic enzymes inside the fruit continued for a certain period of time after harvesting, which may have led to an increase in the levels of these two compounds in dry FLL.Thus, the changes in secoiridoids, in addition to environmental factors, are closely related to their biosynthesis.
The biosynthetic route for the formation of oleosidetype secoiridoids has only been partially reported [23].Based on their structures, it seems that spenuezhenide, oleuropein, nuezhenoside G13, ligstroside, and neonuzhenide were produced from OME-glu esterized with tyrosol, hydroxytyrosol, or their glucosides.OME-glu is a glycosylated product of OME, which is known to mediate oleoside biosynthesis.Glycosylation and esterification have not yet been demonstrated and characterized [4,11].The quantitative analysis of OME and OME-glu showed that there was almost no OME-glu in dry FLL, but a small amount of OME remained.When the content of oleosides such as spenuezhenide and G13 were lower during the six periods, OME-glu also had a lower content, while OME had a higher content (Fig. 1).This variation in the contents of these compounds suggested that the conversion of OME to OME-glu was slow in the later periods (135 DAF to 195 DAF), while the consumption of OME-glu was rapid.During these stages, oleoside content was stable but OME-glu content increased slightly, suggesting that the GT activity of OME on OME-glu decreased gradually in the later harvesting times.To explore the specific GT of OME, we isolated the total FLL protein caused it to react OME.However, the expected chromatographic peak of OME-glu did not appear (Fig. S12).Compared with the deactivation protein, the peak area of OME obviously decreased in the total protein reaction.This may be because the total protein extract was more complex and contained other small-molecule compounds.OME-glu might have reacted rapidly and therefore would not have been detected.According to studies on plant acyltransferases, it is highly probable that enzymes with esterification activity belong to the serine carboxypeptidase acyltransferase (SPCL) family, which utilizes 1-O-β-glucose esters as acyl donors [24].The structure of OME-glu can be recognized as a 1-O-βglucose ester that react with tyrosol, although this has not yet been demonstrated.The esterification activity may also be performed using BAHD, which is named after the first four biochemically characterized ATs, namely benzylalcohol O-acetyltransferase, anthocyanin O-hydroxycinnamoyltransferase, anthranilate N-hydroxycinnamoyl/ benzoyltransferase, and deacetylvindoline 4-O-acetyltransferase [24].Hence, the biological function and biocatalytic potential of esterification and glycosylation warrant further study.These two unknown catalytic steps can potentially be elucidated based on the analysis of the complete genome.In addition, oleoside-type secoiridoids are predominantly found in the Oleaceae family, including plants such as Osmanthus fragrans, olive, and FLL.These plants belong to different genera and exhibit distinct geographical distribution characteristics, yet they have the capability to produce compounds with similar structures.The co-evolution of these plants within the secoiridoids biosynthetic pathway could potentially explain the occurrence of this phenomenon.
Oleoside-type secoiridoids, such as oleuropein and ligstroside, contain phenolic moieties such as tyrosol, hydroxytyrosol, and their C8-glycosylation salidroside and hydroxysalidroside.The production of these compounds is important for the biosynthesis of oleosidetype secoiridoids [15].Tyrosine can react with TyDCs, MAO, and AAS to produce 4-HPAA.Because of the high amino acid sequence identity between TyDC and AAS, they cannot be separated using primary sequence analysis, so them were grouped for analysis [25].There were two expression trends of TyDC/AAS in FLL; therefore, the functions of the sequences require further study.The specific glycosylation step of tyrosol has been explored in R. rosea.Genes homologous to RrT4GT and RrT8GT were cloned and investigated for their potential roles in the glycosylation of tyrosol and hydroxytyrosol (Fig. 5).The maximum likelihood (ML) tree of these homologs was consistent with their protein similarity (Fig. S7D).However, the T4GT homologs, such as Cluster-13715.51696,Cluster-13715.41413, and Cluster-13715.7146, did not have T4GT activities.In contrast, the T8GTs 35,728, 18,726, and 48,377 converted tyrosol to salidroside.35,728, 18,726, and 48,377 were found to be highly expressed in our RNA-seq data (Fig. S7E).Among these, Cluster-13715.18726 has excellent activity, and can completely catalyze tyrosol and hydroxytyrosol.It is unfortunate that the proteins with esterification and glycosylation in the biosynthesis of oleoside-type secoiridoids are still unknown.In order to investigate the biosynthesis of oleosides and characterize the specific enzymes involved, a better understanding of the transition of the compounds between oleosides is essential.This knowledge will enable the reconstitution and optimization of oleoside metabolic pathways, and lays the foundation for the comprehensive development and utilization of FLL to utilize its health-promoting properties.The results of this study provide a basis for reasonable selection of harvest time of FLL, and also provide a foundation for a comprehensive use and development of FLL resources.

Conclusions
In this study, transcriptomic and metabolomic analyses revealed the rules that describe the dynamic changes in the compounds comprising FLL across development periods.Metabolomic data identified 70 compounds across development periods, most of which were secoiridoids.Eleven other compounds were also quantified, such as tyrosol, salidroside, specnuezhenide, oleuropein, G13.The marker compound specnuezhenide was significantly upregulated at 112 DAF, and G13 was enriched at 135 DAF, the coloring period.Oleuropein was abundant at 45 DAF when FLL was in the young fruit stage.Using transcriptome data, we analyzed the biosynthetic pathways of the main components of FLL.In secoiridoid biosynthesis, gene expression was similar between SXS/ OMES and ISY, suggesting that they may be co-expressed that the upregulation of these genes may be beneficial for oleuropein biosynthesis.The in vitro enzyme activity of GTs, similar to T4GT and T8GT, showed their catalytic ability to convert tyrosol and hydroxytyrosol to salidroside and hydroxysalidroside, respectively.Therefore, this study provides valuable information about the harvest period of FLL according to how the compounds accumulate.It also lays foundation for the molecular research on the improvement of FLL quality.

Plant materials and standards
The samples in the present study were harvested from plants at Nanjing University of Chinese Medicine (Nanjing, China), which were identified by the corresponding author of this article (Professor Qinan Wu).Three fruit trees with good growing conditions at the same developmental period were selected, and sampling of fruits was carried out on August 16, 2021 (45 days after flowering, DAF), September,16 2021 (75 days after flowering), October 23, 2021 (112 days after flowering when the color of the exocarp began to change), November 16, 2021 (135 days after flowering), December 16, 2021 (170 days after flowering when the day of the Winter Solstice Festival in China is known to be the best time for harvesting), and January 16, 2022 (195 days after flowering) (Fig. S1).The six sampling times were named as 45 DAF, 75 DAF, 112 DAF, 135 DAF, 170 DAF, 195 DAF, respectively.The harvested samples were quickly frozen in liquid nitrogen and stored at -80℃.The dry FLL served as the standard traditional Chinese medicine which was used as control.Three biological replicates were used for each experiment.

Sample preparation for LC-MS
The biological samples were freeze-dried using a vacuum freeze-dryer (Scientz-100 F).The freeze-dried samples were crushed using a mixer mill (MM 400, Retsch) with zirconia beads for 1.5 min at 30 Hz.Then, 0.3 g lyophilized powder was added to 10 mL 80% methanol solution and treated with ultrasound for 30 min three times.After centrifugation at 12,000 g for 10 min, the extracts were filtered through a 0.22-µm membrane filter before LC-MS/MS analysis.There were six biological repeats in each group.

RNA extension and transcriptome analysis
Total RNA of fruit samples harvested at different times was extracted using the CTAB reagent (QIAGEN, Germany).RNA samples were separately sequenced using an Illumina NovaSeq 6000 platform (Illumina, USA) following the manufacturer's recommendations, in which the mRNA was purified using poly-T oligo-attached magnetic beads.Library quality was evaluated using an Agilent 2100 bioanalyzer.After each library was qualified, the different libraries were pooled according to the

Gene isolation, enzyme expression, and enzymatic activity assay
A HiScript III 1st Strand cDNA Synthesis Kit (+ gDNA wiper) (Vazyme Biotech) was used with 0.5 µg of total RNA to synthesize the first-strand cDNA, which was used as the template for PCR with a 2 × Taq Master Mix (Vazyme Biotech).Gene-specific primers were designed based on annotated results from the transcriptome database.The PCR products were separated using a Gene-JETTM Gel Extraction Kit (Thermo Fisher Scientific), and the purified DNA fragments were cloned into a pET28a vector using a ClonExpress II One Step Cloning Kit (Vazyme Biotech).All clones were sequenced using Sanger sequencing (Tsingke, Beijing, China).The recombinant vector was transformed into Escherichia coli BL21 cells (TIANGEN, Beijing, China) using the heat shock treatment method at 42 ℃ for 45 s.Finally, transformants were plated on LB medium containing 100 ng/L kanamycin.
Single colonies were used to inoculate into LB medium (5 mL) containing 100 ng/L kanamycin.The culture (1 mL) was transferred to 100 mL of the same medium and continued to grow at 37 °C to reach an OD 600 of approximately 1.0.Protein expression was induced by adding isopropyl β-d-thiogalactopyranoside at a final concentration of 1.0 mM.After incubation at 18 ℃ for 20 h, the cells were harvested at 5,000 g, 4 ℃ for 15 min and resuspended in 2 mL lysis buffer (15% v/v glycerinum, 10 mM Tris-HCl [pH 8.0], and 200 mM NaCl).The suspension was then transferred to an ultrasonic cell disruptor system and ultrasonically disrupted for 25 s each, with 35 s on ice between disruptions with three times.The mixture was then centrifuged (30 min, 4 °C, 12,000 g), and the supernatant was carefully collected as the crude protein.
The reaction system of enzyme activity in vitro is as follows: the crude protein (40 µL), 50 mM Tris-HCl (pH 7.5), 500 µM substrate, 2 mM UDPG.They were added to a 250 µL reaction system, and incubated at 30 ℃ in the dark on a shaker.The reaction was stopped by adding 250 µL of methanol.Next, all samples were centrifuged at 12,000 g for 10 min and, the supernatant was subjected to HPLC analysis.

Fig. 1
Fig. 1 The content variation of the main compounds (n = 6).The x-axis shows the six fruit development periods of FLL, with dry FLL representing the traditional Chinese medicine (TCM) of FLL; the y-axis shows the mass percentage of the compound content.Data are represented as mean values with standard deviation obtained from six biological replicates.FLL, Fructus Ligustri Lucidi; DAF, days after flowering

Fig. 3
Fig. 3 Biosynthesis of secoiridoids and expression of related DEGs across the six periods of FLL.Compounds in red were detected in FLL.Red five-pointed stars indicate that the compounds were isomers and could not be distinguished using LC-MS.The words in red indicate compounds that were identified using the standards.The pentagram and the circle-marked compounds indicate that they are isomers.Highly expressed genes are indicated in red, while genes with reduced expression are in blue.The color scales reflect a log 2 -transformed mean of fragments per kilobase of transcript per million mapped (FPKM) values.From left to right, the sample sequence is 45 DAF_1, 45 DAF_2, 45 DAF_3, 75 DAF_1, 75 DAF_2, 75 DAF_3, 112 DAF_1, 112 DAF_2, 112 DAF_3, 135 DAF_1, 135 DAF_2, 135 DAF_3, 170 DAF_1, 170 DAF_2, 170 DAF_3, 195 DAF_1, 195 DAF_2, and 195 DAF_3.FLL, Fructus Ligustri Lucidi; DAF, days after flowering

Table 1
). Tyrosol content increased from 45 DAF to 112 DAF, then decreased from 112 DAF to 195 DAF.Salidroside and hydroxytyrosol accumulated mainly at 45 DAF and 75 DAF, then rapidly decreased at 112 DAF, remaining at stable levels from 112 DAF to 195 DAF.The hydroxysalidroside content increased from 45 DAF to 135 DAF, then decreased gradually from 135 DAF to 195 DAF.The tyrosol, salidroside, hydroxytyrosol, and hydroxysalidroside contents in dry FLL were lower than those at 45 DAF and 75 DAF.Oleoside-11-methyl ester Details of the identified metabolites of FLL unigenes in NR, 24,791 in NT, 8,301 in KO, 5,351 in KEGG, 18,865 in SwissProt, 18,009 a denotes compounds that were compared with the standard substance; Rt represents retention time Table 1 (continued)in PFAM, and 13,371 in GO that were annotated among the DEGs.The number of DEGs were 5,854, 5,467, 1,214, 6,143, and 1,687 between 45 DAF and 75 DAF, 75 DAF and 112 DAF, 112 DAF and 135 DAF, 135 DAF and 170 DAF, and 170 DAF and 195 DAF, respectively.A total of 198 KEGG pathways were mapped, with the most enriched pathways being related to glucose metabolism, fatty acid biosynthesis, and terpenoid biosynthesis (Table